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j ■ Abstract. We report computer simulation results on the cluster formation of dipolar 

Q I colloidal particles driven by a rotating external field in a quasi-two-dimensional setup. 

■ We focus on the interplay between permanent dipolar and hydrodynamic interactions 

. and its influence on the dynamic behavior of the particles. This includes their 

individual as well as their collective motion. To investigate these characteristics, 
we employ Brownian dynamics simulations of a finite system with and without 

' , hydrodynamic interactions. Our results indicate that particularly the translation- 

Q ■ rotation coupling from the hydrodynamic interactions has a profound impact on the 

O ■ clustering behavior. 
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1. Introduction 

Systems of colloidal particles can organize into a wide variety of different structures, 
both under equilibrium conditions ("self-assembly") and in nonequilibrium. Prime 
examples of nonequilibrium structure formation are lane formation of charged colloids 
[HE], shear banding of rod-like particles [3], the coiling up of chains of magnetic colloids 
[1], metachronal waves in driven colloids [5], and the formation of colloidal caterpillars 

Indeed, magnetic particles form a particular subset of colloidal systems, since 
they can be easily manipulated by external magnetic fields. The interplay of particle- 
field interactions and the anisotropic (dipole-dipole) magnetic interactions between the 
particles then leads to a wide variety of structures such as chains [TJ El |9], layers 
P El Uni [m [12], and intricate honeycomb-like structures p^lTil fTH]. 

In this study we want to focus on the pattern formation of dipolar particles that are 
driven by a rotating field. In three-dimensional systems dipolar particles exposed to such 
a field tend to form layers (biaxial field) [161 E] or membrane-like structures (triaxial 
field) [13l[14]. In a two-dimensional geometry, in which the external field rotates in the 
plane of the particles, the dipoles tend to agglomerate into clusters. This phenomenon 
has been observed experimentally [T71 |T8] as well as in computer simulations that were 
performed by two of the authors of this study ^9\. Clustering does not only occur for 
particles carrying a permanent dipole moment, but also for particles with an induced 
dipole moment [2011211122]. 

At suitable fields strengths and frequencies, the (permanently) dipolar particles 
perform a synchronous rotation with the field. This motion gives rise to an attractive, 
long-range particle interaction (oc — 1/r^), which induces a first-order phase transition 
between a dilute and a denser phase (in an infinitely extended system). The observed 
cluster formation then corresponds to spinodal decomposition inside of the coexistence 
region [19] . 

It has been shown that hydrodynamic interactions play a crucial role in numerous 
colloidal systems [23] |2H [25] [26] . The influence of hydrodynamic interactions on the 
clustering phenomenon was, however, only shortly touched upon in [19]. Here, in this 
study, we want to look into the effects that an implicit solvent has on the cluster 
formation and the dynamic behavior of the particles in more detail. 

The main tool we use are Brownian dynamics computer simulations with and 
without hydrodynamic interactions. The particles are modeled by a quasi-two- 
dimensional system of soft spheres with permanent dipole moments. This means that 
the dipoles can rotate freely in all the spatial directions, while the translational motion 
is restricted to a two-dimensional plane. In our hydrodynamic simulations, we not only 
account for the translational hydrodynamic couplings, but for all the couplings between 
the translational and rotational motions of the particles. To understand the influence 
of these different couplings is one of the goals of this study. Contrary to [19], we here 
consider a finite system size, i.e., we do not employ periodic boundary conditions. 
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This paper is organized as follows: After introducing the model and the simulation 
techniques, we discuss the influence of hydrodynamic interactions on the regions of 
cluster formation in the field strength-frequency domain. To understand the differences 
stemming from the solvent interactions, we also consider the magnetization and 
synchronization behavior of the particles. In a next step, we investigate the dynamics 
of the clusters and their formation. Finally, we examine the structure of the clusters 
with respect to their internal order. The paper is then closed with a brief summary and 
conclusions. 



2. Model and simulation methods 



In this study, we consider a quasi-two-dimensional system of dipolar colloidal particles 
that are immersed in a solvent. As a model for the colloids we use a dipolar soft sphere 
(DSS) potential, which is comprised of a repulsive part U'^^^ and a point dipole-dipole 
interaction part U^: 

t/°ss(r,„ M., h) = U'^^r,,) + t/°(r,„ /x,, fi,) (1) 

In ([T]), Tij is the vector between the positions of the particles i and j, rij its absolute 
value, and is the dipole moment of the zth particle. The dipolar and repulsive 
interaction potentials are given by 



K'-iji P-i, H'j) — 5 3 V^J 

ij ij 

and [271 

dr/ss 

f/-P(r) = [/SS(r) - f/SS(r,) + (r^ - r)^(r,), (3) 
respectively. Here, U^"^^ is the shifted soft sphere potential, where 

U^^[r)=Ae(—Y (4) 



is the unshifted soft sphere (SS) potential for particles of diameter a. 

We investigate the system by making use of Brownian dynamics (BD) simulations 
with and without hydrodynamic interactions (His). Specifically, the positions of the 
particles are evolved in time via [281 l29l [30] 



,(t+/\t) = v,{t) + ^Y.^j^^,M 

+ rTYl D^T, At + R,(D, At) (5) 
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whereas their orientations = Hi/fi evolve according to 

e,(t + At) = e,(t) + 

+ YfYl ^fj^At) X e,(t) + R,(D, At) X e,(t). (6) 

The forces and torques in ([5]) and ([6]) are given by 

F, = -Vr,$^f/°^'(r,„M.,M,), (7) 
Ti = T°ss + Tf*, (8) 

where 

s = X V^, M., M,), (9) 

= tiiX B"^*. (10) 

In ( fTOl) . B'^^* denotes an external field that is homogeneous and rotates in the plane of 
the dipolar monolayer. Specifically, 

B""""* (t ) = i?o ( cos + sin Wot ) , (11) 

where Uq is the frequency of the field and Bq its strength. 

In (El) and (El), D^t, d^.R, D^^^^ ^^^^ j^rr subtensors of D^^, D™, D^^, and 
are subtensors of the grand diffusion tensor 

and will be specified below. The random displacements in (jSj) and (jSj) behave according 
to 

(Ri) = 0, (13) 
(Ri(At)Rj(At)) = 2'DijAt. (14) 

The actual calculation of these displacements can be done by evaluating 

R=y2AtL-^, (15) 

where R is the vector comprised of all the Rj, ^ is a vector of normally distributed 
random numbers, and L is a lower triangular matrix which satisfies 

D = L ■ L^. (16) 
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In the present study, we take the His into account up to third order in the inverse 
particle distance, which corresponds to a far field approximation. The tensors D"""""", 
j^TR^ qrt^ g^^^ qrr g^^g ^i^^^ gi^gj^ bv [29| [301126] 



rprp 



D,^ = -p^I, (17) 



TT 

ksT 1 



2o-2 

(I "I" ^ij^ij) "I" o 2 3f jj,Tjj^ 



mrj rij 

= D™ = 0, (19) 

D™ = D^^t ^ |£^^ef„ (20) 
Snrj r^j 

nr^s^i, (21) 



where is the viscosity of the solvent and e the Levi-Civita density, which satifies 
ei23 = ^231 = ^312 = 1, £321 = 6213 = €132 = " 1, and is equal to zero for other choices 
of the indices. These tensors describe the different hydrodynamic couplings, i.e., the 
coupling between the translational motion of the particles (TT), between the translation 
and the rotational motion (TR) and vice versa (RT), and the coupling between the 
rotational motion (RR). The specific tensor given in dlTl) and ([18]), i.e., D^^, is the well 
known Rotne-Prager tensor [21]. 

Note that by setting to zero all the hydrodynamic coupling tensors involving 
different particles i j in ([T711 - ([22ll . one arrives at the standard algorithm for BD 
simulations without His [32]. 

Considering the equation of motion ([5]), we can see that the TR coupling relates all 
the torques acting on the particles, i.e., their rotational motions, to all the translational 
motions of the particles. Physically, the TR coupling describes how the translational 
motion of the particles is affected by the flow flelds in the solvent that are caused by 
the rotations of the particles. Indeed, as we will see later, the TR/RT coupling is 
particularly relevant for the overall dynamical behavior of the system. It is therefore 
instructive to briefly illustrate the implications of this coupling for a simple two-particle 
system (for details, see [30l [26] ). 

To this end, let us consider two particles that are located at a distance from each 
other on the x-axis in a right-handed coordinate system. Due to the TR coupling, 
anticlockwise rotation (following the application of a torque) of the particle at the 
larger value of x results in the other particle moving in the negative y-direction. By 
realizing that the flow flelds follow the rotation of the particle, this process can be easily 
understood. On the other hand, an anticlockwise rotation of the particle at smaller 
values of x causes the other particle to move into the positive y-direction (for a sketch, 
see [30]). 

In this study, we consider = 324 particles in a simulation box that is bounded 
by soft walls (cf. |33]), i.e., we do not use periodic boundary conditions. Therefore, it 
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is not necessary to use special techniques (e.g., Ewald sums) to treat the long-range 
dipole-dipole interactions. The forces and torques can be calculated directly via ([2]). 

For convenience, we make use of the following reduced units: Field strength 
Bq = {a^/eY^^Bo; frequencies of the field Uq = Uoa'^/D'^, where = kBT/Sirrja is 
the translational diffusion constant; dipole moment /i* = (ecT^)~^/^/i; time t* = tD^ /a'^; 
temperature T* = ksT/e; position r* = r/a. In the following we specialize to systems 
at temperature T* = 1 and of dipole moment fi* = 3. This choice corresponds to a 
dipolar coupling strength of A = /i*^/T* = 9 that is sufficiently large to enable the 
system to form clusters for suitable field strengths and frequencies [19]. The density of 
the particles in the simulation box is of no importance in the investigated systems, since 
the particles typically agglomerate into a single cluster. 

3. Dynamics on the particle level 

Applying a rotating field in the plane of a dipolar monolayer can cause the dipolar 
particles to agglomerate into two-dimensional clusters [I9] . This agglomeration is caused 
by a synchronization phenomenon: At suitable field strengths and frequencies, the 
particles follow the field at the same phase difference resulting in an effective interparticle 
interaction of the form [SU [19] 

U'^'ir,) = (23) 

Equation (!23|) can be arrived at by averaging the dipole-dipole potential over one 
rotational period of the field under the assumption that /Xj(t) = l-ijit) = fio[Gx cos{uot + 
6) + eySm{uJot + 6)]. Here, 6 is some phase difference. A further crucial assumption 
in the derivation of ( 123|) is that the translational motion of the particles during one 
rotational period of the field is negligible. 

The potential fl23|) is attractive, which leads to the aforementioned cluster 
formation. Moreover, in an infinitely extended system, i.e., in the thermodynamic limit, 
the potential f l23|) gives rise to a first order phase transition [T9] . 




Figure 1. Snapshots of systems at = 50 and loq ~ 240 after clusters have formed, 
(a) Without and (b) with His. 
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Simulation snapshots of systems after clusters have formed can be seen in figures 
[U^a) and (b). The former shows the snapshot of a system [B^ = 50, Uq = 240), whose 
particles do not interact via His, while His are included in the system associated with 
the latter snapshot. Since the formation of clusters can be observed in both the systems, 
we note as a first result that cluster formation is not prevented by the presence of His. 
As already noted in [19], this is not a priori clear, since His can induce additional 
motion in a nonequilibrium system and the averaged potential ( 123|) is only valid as an 
approximation to the true interparticle interaction if the translational motion of the 
particles during one rotational period of the field is small. This is certainly the case for 
the system without His: Inspecting the mean squared displacements (averaged over all 
the particles) at Bq = 50 and Uq = 240, we find that the field rotates about 30 times 
before a particle transverses a distance of its own diameter a. 

The driving frequency Wq = 240 chosen in figure [U corresponds to Uo ~ 54 MHz, 
if we assume the values of the diffusion constant (Dq ~ 38 /xm^/s) and particle size 
{a = 13 nm) that are given in [29] for a ferrofluid. 

Note that the cluster formation in infinitely extended quasi-two-dimensional 
systems corresponds to spinodal decomposition within the coexistence region of a phase 
transition [19]. This is not exactly true for the cluster formation we observe in the 
present study. The finite systems considered here do not undergo a phase transition. 
However, the used system is well suited as a model system to investigate the influence 
of the His on the dynamic behavior of the individual clusters. 

In the following we ask how the collective rotational behavior of the particles 
changes if solvent-mediated interactions are taken into account. First, consider figure 
[21 which shows whether cluster formation occurs for selected state points in the field 
strength-frequency domain. Presented are results for both the cases with and without 
His included. Compared to the simple BD system cluster formation breaks down at 
smaller frequencies Uq in the hydrodynamically interacting system. In the former, cluster 
formation can be observed up to Uq ~ 450 (at -Bq = 50) while cluster formation ceases 
at Uq ~ 350 when His are present. 

To understand the breakdown of cluster formation in more detail, we now examine 
the rotational motion of the particles. As explained above, synchronized rotation is 
necessary for (l23i) to hold, i.e., is a prerequisite for cluster formation. Figure |3] shows 
the absolute value of the magnetization normalized with respect to its saturation value 



over the driving frequency uJq of a system [Bq = 50, cf. Fig. [2] with all the hydrodynamic 
couplings and without His included. The magnetization indicates how aligned the 
particles are in a given state, i.e., indicates if they follow the field. As can be seen, 
the magnetization starts with values close to one for both the systems, corresponding 
to an aligned state. At Uq ~ 270 the magnetization begins to drop for the system that 




(24) 
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Figure 2. Cluster formation (top) with and (bottom) without His in the field strength- 
frequency domain. Squares/circles indicate where cluster formation occurs/does not 
occur. The systems are at temperature T* = 1 and dipole moment /x* = 3. 



includes His. The magnetization in the system without His remains at M/Mq 1 up 
to Lu^^ ^ 420. 
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Figure 3. Magnetization normalized with respect to its saturation value over the 
driving frequency of the external field for a system with and without His. The fields 
used are of strength Bq = 50. The vertical line indicate where cluster formation ceases 
with (dashed) and without (solid) His. 



This magnetization behavior implies that the particles in the hydrodynamically 
interacting system are less ahgned with each other for co^ > 270. In particular, the 
synchronization breaks down at lower frequencies resulting in a premature breakdown 
of cluster formation. 
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4. Cluster dynamics 

We now turn to the dynamics of the entire cluster. In [19j , it was already shown that in 
systems with His, clusters form at a significantly faster pace (as compared to systems 
without His) . In this section, we aim to investigate the changes induced by His in more 
detail. 
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Figure 4. Mean distance between tlie particles over time d in a system {Bq = 50, uJq = 
240) with all hydrodynamic interactions included (All HI), without the hydrodynamic 
RT/TR coupling (No RT/TR), without the hydrodynamic RR couphng (No RR), and 
without aU His (No HI). 

First, consider figure HJ which shows the mean distance 

9 ^ 

^>) = Miv3T)5:E^«.W (25) 

between the particles for a system at Bq = 50 and Uq = 240. These particular field 
parameters were chosen for three reasons: First, the frequency is sufficiently high 
to ensure that the effective potential fl23|) describes the interparticle interaction well. 
Second, as seen in figure [21 cluster formation occurs for both a hydrodynamically as 
well as a not hydrodynamically interacting system. Third, the magnetization of the 
systems with these choices of Uq and Bq is maximal and identical irrespective of the 
presence of His [cf . figure [3] . 

Specifically, figure H] shows the evolution of d over time for particles interacting via 
His including all the hydrodynamic couplings, for particles lacking the hydrodynamic 
RT/TR coupling, particles lacking the RR coupling, and particles not interacting via 
His at all. In all the cases, d assumes a constant minimal value at long times. To 
understand this, recall that we consider a single simulation box filled with particles here. 
In an infinite system, the cluster would keep growing in time indefinitely with a power 
law behavior [351 ESI EZ] , since the process corresponds to spinodal decomposition [19j . 
Here, however, the growth process stops once all the particles have been incorporated 
into the cluster and a stationary state is reached. 
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In the systems that include solvent-mediated interactions (All HI, No RT/TR, 
No RR), the value of d drops significantly faster than in the case without any His. 
Consequently, the average distance between the particles decreases faster, which means 
that the cluster formation process is sped up. The acceleration is neither infiuenced by 
the lack of the RT/TR nor the RR coupling, which implies that the TT coupling alone 
is responsible for this effect. Note, however, that the lack of the presence of the RT/TR 
coupling expresses itself by a different value of d at long times [see figure H] . 
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Figure 5. Mean angular frequency uj*^ around the cluster center of the particles over 
distance from particle center r* for a system at Bq = 50 and uJq = 240. Shown are 
values for a system interacting via all the hydrodynamic couplings and for a system 
lacking the RT/TR coupling. 



The RT/TR coupling does have another interesting influence on the dynamic 
behavior of the particles. In figure [5l the mean angular frequency of the particles with 
respect to the cluster center over the distance from the center is shown for a system 
at Bq = 50 and Uq = 240. Values for the hydrodynamically interacting case with all 
the couplings included as well as the case lacking the RT/TR coupling are presented. 
In the system that includes the RT/TR coupling, the angular velocity of the particles 
differs from zero at all the displayed distances from the cluster center r*. Hence, the 
particles perform a rotation around the cluster center in the rotational direction of the 
external field. The system lacking the RT/TR coupling does not show such a rotational 
behavior. As can be seen in figure Ej the mean angular frequency of the particle around 
the cluster center is essentially zero at all distances. 

This collective rotation is caused by the individual, field-driven rotations of the 
particles. The rotational motion of the particles creates a fiow field that induces 
translational motion in all the other particles [see the argument given below f|T7|) - fl22|) ]. 
Therefore, the TR coupling alone is responsible for this behavior. The RT coupling does 
not contribute in any way to the cluster rotation. 

As we have seen in the previous section. His result in a premature breakdown of 
cluster formation, i.e., a breakdown at smaller driving frequencies of the field (relative 
to the case without His). It stands to reason that the cluster rotation induced by 
the His has a significant infiuence on this behavior. The particles perform additional 
translational motion (around the cluster), which makes the effective potential ( 123|) less 
accurate as a description for the interparticle interaction at a given driving frequency. 
The more the particles move during one period of the field, the less does the effective 
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potential capture the actual interaction between the particles. 

Finally, note that the RR coupling does not seem to have any significant infiuence 
on the dynamic behavior of the cluster. That is, it does not contribute to the accelerated 
cluster formation or the cluster rotation. The former is illustrated by figure HI which 
shows that the mean distance d behaves essentially identically to the system with all the 
His included. The fact that the cluster rotation is not influenced by the RR coupling 
can be seen in flgure [51 Despite the presence of the RR coupling, the cluster does not 
rotate if the RT/TR couphng is absent. 

5. Internal structure of the cluster 

In a recent experimental study, Weddemann et al. [IT] showed the existence of cluster 
formation in two-dimensional systems of (permanently) dipolar particles that are driven 
by a rotating external fleld. In particular, the authors of [IT] observed the formation of 
hexagonally ordered particle agglomerates in their experiments. 

The general clustering phenomenon of dipolar particles exposed to rotating fields 
was later identified as spinodal decomposition in a simulation study by two of the 
authors of the present study [19]. However, to reproduce clusters of hexagonal order, 
very low temperatures inside of the two-phase coexistence region had to be considered 
in computer simulations. It was conjectured that the hexagonally ordered clusters occur 
in the vapor-solid coexistence region, i.e., at coupling strengths above the ones related 
to the vapor-liquid region. 

As shown in the previous section. His can have a significant infiuence on the 
collective motion of the particles. Here, in this section, we want to investigate, whether 
His preserve the internal cluster structure. Despite experimental evidence of the 
hexagonal order, this fact is debatable since it remains unclear to what degree the 
rotational motion of the magnetic particles in the experimental work [17] follows the 
dipole moment. 

In order to gain insight into the emergent (hexagonal) structures in the present, 
finite systems, we consider the bond order parameter 



at different dipolar coupling strengths A. Here, Afn is the set of nearest neighbors of 
particle n, which consists of particles that are closer to particle n than the distance of 
the first minimum in the pair correlation function of the system. The systems considered 
in the following are of sufficiently high coupling strength ensuring that cluster formation 
does indeed occur. 

In figure |6l ipQ as function of the dipolar coupling strength A for systems {Bq = 50) 
with and without His at oJq = 240 and 300 is shown. We note that ipQ is essentially 
independent of time in the stationary situation, where all the particles are part of 
the cluster. At cJq = 300, the bond order parameter increases with A for both the 




(26) 



Dynamics of cluster formation in driven dipolar colloids dispersed on a monolayer 12 



to 



l.Or 

0.8 
0.6 
0.4 

0.0^ 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



2 5 



= 240 



2 2 



n O ° ° 
n O O 

2 ° A * ^ - * 


O o 

A A A 


Wq* =300 


O O No HI 
A A All HI 



10 



12 



14 



16 



Figure 6. Bond order parameter ipQ over coupling strength A for a system {Bq — 50, 
T* = 1) with and without His at (top) loq — 240 and (bottom) cjq = 300. 



hydrodynamically as well as the not hydrodynamically interacting system. If His are 
not present, an increase from ipQ ^ 0.52 at A = 7.02 to ipQ ~ 0.72 at A = 16 can 
be observed. This increase qualitatively agrees with the one observed in the Langevin 
dynamics simulations in [19] . In the system with His included, the bond order parameter 
V'e increases considerably less with increasing coupling strength. However, there is still 
significant order in the system. 

In the system with Uq = 240 and His included, ipQ behaves similarly to the case 
with Uq = 300 and His. If these interactions are not taken into account, however, less 
order than in the Uq = 300 system can be observed. The smaller frequency allows for 
more translational motion during one rotational period of the field, resulting in more 
spatial inhomogeneity. 

Now consider figure [71 which shows the pair correlation functions of a system 
{Bq = 50, Uq = 300) with and without His included. The pair correlation function of 
the not hydrodynamically interacting system shows a double peaked maximum after the 
first minimum, which is typical for a hexagonally ordered system. The hydrodynamically 
interacting system, on the other hand, does not have this feature. Hence, the particles 
tend to have six angularly equally distributed neighbors as shown by the value of ipe 
[cf. figure [6], but seem to lack the long-range positional order of a hexagonally structured 
system. Further, the extrema are much more pronounced in the system without His, 
indicating a more ordered state. 

In conclusion, in both systems in figure O the His tend to weaken the hexagonal (or 
hexatic) order present in the system. The lower value of ipQ in systems with His can be 
explained by the collective cluster rotation induced by the hydrodynamic TR coupling. 
The particles rotate around the cluster center and do not stay at fixed lattice sites. This 
behavior results in a reduction of the bond order parameter and the hexagonal structure 
in the system. 
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Figure 7. Pair correlation functions of systems at Bq — 50 and loq — 300 with and 
witliout His included. 



6. Conclusions 

In this study, we liave investigated the influence of the solvent on the dynamic behavior 
of rotationally driven dipolar particles. 

The cluster dynamics and formation is influenced in three major ways by His. First, 
His accelerate the cluster formation process. While the particle approach each other, 
flow flelds are created that drag other particles along. Second, we have shown His to 
induce a collective cluster rotation. Without His no such rotation can be observed. The 
driven rotating particles pull the solvent with them in their rotational motion, which 
results in the translational motion of the particles around the cluster center. As a last 
major point, we have shown cluster formation to cease at lower driving frequencies of 
the field. We attribute this to an earlier breakdown of synchronization if His are present 
and an increase in translational motion due to the collective cluster rotation. 

Moreover, we have studied the influence of His on the internal structure of the 
clusters. Our results indicate that His tend to decrease the hexagonal order in the 
systems. 

Recently, one of the authors of this study has conducted an investigation of the 
structure of a closely related system [38] . The system consisted of particles that only 
interact via His, are conflned to a monolayer, and have flxed angular velocities. It was 
found that hexagonal particle agglomerations rotate, with the hexagonal order melting 
and recrystallizing periodically. In the system focused on in the present study, no such 
phenomenon could be observed. We attribute this to the presence of dipolar interactions 
in our system and a lack of Brownian motion in [3H] • 

Another study has recently investigated the interplay of conflnements and His |39] . 
It was shown that particle rotations can be utilized to create directed translational 
motion in colloids in speciflc geometries. This result suggests that it might be interesting 
to examine the effects of different conflnements on the dipolar systems considered here. 

In conclusion, we have shown that His have a considerable influence on the 
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formation and the dynamics of clusters of driven dipolar colloidal particles in a two- 
dimensional geometry. Usually, His seem to affect colloidal nonequilibrium clustering 
phenomena less significantly than in the present system. As an example, consider the 
process of colloidal gelation in two-dimensional Lennard- Jones systems [IQ]. In contrast 
to our system, the agglomeration of the particles in [40] is only marginally influenced by 
the His. In general, however. His can significantly alter nonequilibrium processes. For 
instance. His can enhance ratchet effects |¥T| H2] or synchronize the motion of eukaryotic 
[13] or bacterial jB] fiagella. 
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